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Abstract. Critical behaviour of a nearly critical system, subjected to vivid turbulent 
mixing, is studied by means of the field theoretic renormalization group. Namely, 
relaxational stochastic dynamics of a non-conserved order parameter of the Ashkin- 
Teller-Potts model, coupled to a random velocity field with prescribed statistics, is 
considered. The mixing is modelled by Kraichnan's rapid-change ensemble: time- 
decorrelated Gaussian velocity field with the power-like spectrum oc k~ d ~^. It is shown 
that, depending on the symmetry group of the underlying Potts model, the degree of 
compressibility and the relation between the exponent £ and the space dimension 
d, the system exhibits various types of infrared (long-time, large-scale) scaling 
behaviour, associated with four different infrared attractors of the renormalization 
group equations. In addition to known asymptotic regimes (equilibrium dynamics of 
the Potts model and the passively advected scalar field), existence of a new, strongly 
non-equilibrium type of critical behaviour is established. That "full-scale" regime 
corresponds to the novel type of critical behaviour (universality class) , where the self- 
interaction of the order parameter and the turbulent mixing are equally important. 
The corresponding critical dimensions depend on d, £, the symmetry group and the 
degree of compressibility. The dimensions and the regions of stability for all the regimes 
are calculated in the leading order of the double expansion in two parameters £ and 
e = 6 — d. Special attention is paid to the effects of compressibility of the fluid, because 
they lead to nontrivial qualitative crossover phenomena. 



PACS numbers: 05.10.Cc, 64.60.Ht, 64.60. ae, 47.27.ef 
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1. Introduction 

Numerous complex systems, involving "infinitely" many degrees of freedom, 
demonstrate very interesting behaviour in the vicinity of their critical points. The 
typical feature of such systems is the divergence of the correlation length, although the 
underlying microscopic interactions are short-range ones. As a result, the behaviour 
of such systems "decouples" from their specific physical nature (fluids, magnets, 
superconductors etc.) and becomes highly universal: their thermodynamic and 
correlation functions exhibit scaling (self-similar) behaviour and depend only on "global" 
characteristics of the system (like its symmetry and dimension). Thus one can speak 
about the theory of critical state (or critical behaviour), in contradistinction with the 
numerous special theories of diverse critical phenomena (of course, very interesting and 
important on their own). 

The adequate theory of critical state (both qualitative and quantitative) is based 
on the field theoretic renormalization group (RG); see the monographs [U [2] for the 
detailed presentation and the bibliography. In the RG approach, possible types of critical 
behaviour (universality classes) are associated with infrared (IR) attractive fixed points 
of certain renormalizable field theoretic models. 

Most typical equilibrium phase transitions belong to the universality class of the 
0(n)-symmetric if) 4 model of an n-component scalar order parameter if). Universal 
characteristics of the critical behaviour (critical exponents, amplitude ratios, scaling 
functions) depend only on n and on the spatial dimension d and can be calculated within 
systematic perturbative schemes, in particular, in the form of expansions in e = 4 — d, 
the deviation of the space dimension from its upper critical value d = 4. 

More exotic, but still actual and important example is provided by the Ashkin- 
Teller-Potts (or simply Potts) class of models [3] [9]. In rough words, such a model 
describes a certain system which locally has n states, while the energy of any given global 
configuration depends on whether the pairs of neighboring sites are in the same state or 
not. In the continuous formulation, Potts-type models are conveniently represented 
by the effective Hamiltonian for the n-component order parameter with a trilinear 
interaction term, invariant under the n- dimensional hypertetrahedron symmetry group 

mm- 

The Potts-type models have numerous physical applications: solids and magnetic 
materials with nontrivial symmetry, spin glasses, nematic-to-isotropic transitions in the 
liquid crystals, bond percolation problem, random resistor networks and many others, 
see the papers [5] [S] and the references therein. The Potts model (especially its 
paradigmatic special case n = d = 2, the famous two-dimensional Ising model) has 
long become a source of inspiration for the new physical and mathematical ideas like 
integrability, conformal invariance, discrete holomorphisity etc.; for a recent discussion, 
see the papers [10], materials of the meeting [11] and references therein. 

The problem of the nature of the phase transition in the Potts model has a long and 
entangled history; see e.g. the discussion in [Sj. According to naively applied Landau's 
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phenomenological approach, existence of a triple term excludes the possibility of the 
second-order transition. On the contrary, exact two-dimensional results, numerical 
simulations and RG analysis suggest that for n small enough, the phase transition 
in the Potts models is of the second order. We will not try to shed new light onto 
this interesting but difficult question: it should first be solved for the original static 
problem. To proceed, in this paper we accept the point of view that the existence of 
an IR attractive fixed point of the RG equations implies existence of a certain scaling 
(self-similar) IR asymptotic regime and, therefore, the existence of some kind of critical 
state. 

It has long been realized that the behaviour of a real system near its critical point 
is extremely sensitive to external disturbances, gravity, geometry of the experimental 
setup, inclusion of impurities and so on; see the nice book [12] for the general discussion 
and references. The ideal critical behaviour of a pure equilibrium infinite system can 
be distorted by finite-size effects, limited accuracy of measuring the temperature, finite 
time of evolution (ageing) etc. What is more, some disturbances (randomly distributed 
impurities in magnets and turbulent mixing of fluid systems) can change the very type 
of the phase transition (first-order to second order one, and vice versa) and produce 
novel types of critical behaviour (universality classes) with rich and rather unexpected 
properties. 

Investigation of the effects of various kinds of deterministic or chaotic flows on 
the behaviour of the critical fluids (like liquid crystals or binary mixtures) has shown 
that the flow can destroy the usual critical behaviour: it can change to the mean-field 
behaviour or, under some conditions, to a more complex behaviour described by new 
non-equilibrium universality classes |13j-[22j. 

In this paper we apply the field theoretic RG to study the effects of turbulent mixing 
on the critical behaviour of the (generalized) Potts model. Special attention will be paid 
to compressibility of the fluid, because it can lead to interesting crossover phenomena. 

Bearing in mind applications to liquid crystals or percolation in moving media, we 
consider a purely relaxational stochastic dynamics of a non-conserved order parameter, 
coupled to a random velocity field. 

The latter is modelled by the Kraichnan's rapid-change ensemble: time-decorrelated 
Gaussian field with the pair correlation function of the form (vv) cx S(t — t')k~ d ~^, 
where k is the wave number and < £ < 2 is a free parameter with the most realistic 
( "Kolmogorov" ) value £ = 4/3. Some time ago, the models involving passive scalar 
fields advected by such "synthetic" velocity ensembles attracted a great deal of attention 
because of the insight they offer into the origin of intermittency and anomalous scaling 
in the real fluid turbulence; see the review paper [23] and references therein. The RG 
approach to the problem of passive advection is reviewed in |24j. 

In the context of our study it is especially important that Kraichnan's ensemble 
allows one to easily model compressibility of the fluid, which appears rather difficult if 
the velocity is modelled by dynamical equations; cf. [21]. 

The outline of the paper is the following. In section [2] we give the detailed 



Critical behaviour of the stirred dynamical Potts model 



4 



description of the model and its field theoretic formulation. Following [9] , we consider the 
generalized model with a certain symmetry group Q; the plain Potts model corresponds 
to the group Z n . In section [3] we analyze the ultraviolet (UV) divergences and show 
that the model is multiplicatively renormalizable. Thus the RG equations are derived 
in the standard fashion; see section HI The corresponding RG functions (/3 functions 
and anomalous dimensions 7) are calculated in the one-loop approximation. 

In section [5] we show that the model reveals four different types of critical behaviour, 
associated with four possible fixed points (strictly speaking, attractors of a more 
general form) of the corresponding RG equations. Three regimes correspond to known 
situations: Gaussian or free field theory, passively advected scalar field without self- 
interaction (the nonlinearity in the original stochastic equation is irrelevant), and the 
original Potts model without mixing. The most interesting fourth regime corresponds 
to a novel non-equilibrium universality class, formed by the interplay between the self- 
interaction and turbulent mixing. 

A given critical regime can realize if the corresponding fixed point is admissible 
from the physics viewpoints: it should be IR attractive and lie in the physical range 
of model parameters. In our problem, admissibility of the fixed points depends on the 
spatial dimension d, exponent £, symmetry group Q and the degree of compressibility. 
The regions of admissibility form a pattern in the space of model parameters, much 
more complicated in comparison with the ?/> 4 model in the same velocity ensemble [21] 
or the Potts model in a strongly anisotropic shear flow [22]. As a result, interesting 
crossover phenomena occur as the degree of compressibility is varied. These issues are 
discussed in great detail in the same section |5j 

In general, the critical dimensions in our problem are calculated as double series in 
e = 6 — d and £. In section [6] we present the explicit first-order (one-loop) expressions 
for the critical dimensions of the basic fields and parameters. Section [7] is reserved for 
a brief conclusion. 

2. Description of the model. Field theoretic formulation 

Relaxational dynamics of a non-conserved n-component order parameter ip a (x) with 
x = {t, x} is described by the Langevin-type stochastic differential equation 



where d t = d/dt, Ao > is the (constant) kinetic coefficient and ( a (x) is a Gaussian 
random noise with zero average and the pair correlation function 



where d is the dimension of the x space. Here and below, the bare (unrenormalized) 
parameters are marked by the subscript "0." Their renormalized counterparts (without 
the subscript) will appear later on. The normalization in (12. 2p is chosen such that the 




(2.1) 




(2.2) 
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steady-state equal-time correlation functions of the stochastic problem are given by the 
Boltzmann weight exp{— T-L(ip)}. 

Near the critical point, the static Hamiltonian 1-L{ip) of the Potts model is taken in 
the form El d 

W)= j &Ui(x)9 2 i(x) + |i(x)i(x)+ 

+ !yiU c Va(x)V&(x)V>c(x)}, ( 2 ' 3 ) 

where di = d/dxi is the spatial derivative, d 2 = didi is the Laplace operator, 
t oc (T — T c ) measures deviation of the temperature (or its analog) from the critical 
value and go is the coupling constant. Summations over repeated indices are always 
implied (a,b,c = 1, . . . , n and i = 1, . . . , d). After taking the functional derivative 
6T-L(ip) /"0(x) one has to replace r/)(x) in (12. ip by the time- dependent field ip{x). 

Following |9j, we consider the generalized case of certain symmetry group Q, which 
has the only irreducible invariant third-rank tensor R a b c ', without loss of generality it 
is assumed to be symmetric. In the one-loop approximation we will only need to know 
the coefficients in the contractions 

RabcRabe = Rl^ce, RaecRchbRbfa = RlRehf- (2-4) 

In the original Potts model Q = Z n , the symmetry group of the hypertetrahedron 
in the n-dimensional space. Then it is convenient to express the tensor R a b c in terms of 
the set of (n + 1) vectors e a which define its vertices Ej: 

F? l — V 1 p a p a p a 

1 L abc / j °a °fe °c i 



where the e° satisfy 

n+1 n+1 n 

a=l a=l a=l 

For the coefficients in (12.4=1) the relations (12.51) give 

R 1 = (n+ l) 2 (n- 1), i? 2 = (n + l) 2 (n-2). (2.6) 

The stochastic problem (12. ip . (12. 2p can be reformulated as the field theoretic model 
of the doubled set of fields $ = with the action functional 

<S(V, V f ) = AoVM+V'l (~3t + A 9 2 -A r ) iPa-goXoR^J^. (2.7) 

Here ^' = ^'(t, x) is an auxiliary "response field" and all the needed integrations over 
the arguments of the fields and summations over the repeated indices are implied, for 
example 



ipld t ip a = ^2 \ dt \ dxipl(t,x)d t ip a (t,x). 

a=l J J 



This formulation means that the statistical averages of the random quantities in 
the original stochastic problem can be represented by functional integrals over the full 



Critical behaviour of the stirred dynamical Potts model 



6 



set of fields with the weight exp«S($), and can therefore be interpreted as the Green 
functions of the field theoretic model with the action ( 12. 7p . 

The model ( 12. 7ft corresponds to the standard Feynman diagrammatic technique 
with two bare propagators (if)if>>)o, (if>if>)o and the trilinear vertex ~ ip^ip 2 . In the 
frequency-momentum (u;-k) representation the propagators have the forms 

<V.lfc>o(«, k) = - 2 2 Ao ff 2 , (2.8) 
u 2 + Ag (k 2 + r ) 

where k — |k| is the wave number. 

The Galilean invariant coupling with the velocity field v = {vi(t, x)} for the 
compressible fluid (c^ ^ 0) is introduced by the replacement 

d t ip -»■ d t 4> + a + (ao - l)(vidi)if) = Vi^ + a (<9^)^ (2.9) 

in (12.11) . Here Vt = dt + UjcJj is the Galilean covariant (Lagrangian) derivative and 
ao is an arbitrary parameter. For the linear advection-diffusion equation, the choice 
ao = 1 corresponds to the conserved quantity if) {e.g. density of an impurity) while 
ao = is referred to as "tracer" (concentration of the impurity or the temperature of 
the fluid; in this case, the conserved quantity is if)'). In the presence of a nonlinearity 
in (12. ip . it is necessary to keep all the terms in (I2.9P in order to ensure multiplicative 
renormalizability |21j . 

In the real problem, the field v(t, x) is governed by the Navier-Stokes equation. 
Here we employ a simplified rapid-change model, in which the velocity obeys a Gaussian 
statistics with zero average and the prescribed correlation function 

(vi(t, x)vj{t', x')> = 5(t-t') Aj(r), r = x-x' (2.10) 

with 

(27p ¥H { P «M + «^( k )> exp(ikr). (2.11) 

Here -Pij(k) = 5^- — kikj/k 2 and Qij(k) = kikj/k 2 are the transverse and the longitudinal 
projectors, respectively, D > is an amplitude factor and a > is an arbitrary 
parameter which measures the degree of compressibility. The case a = corresponds 
to the incompressible fluid (diVi = 0), while the limit a — > oo at fixed aD corresponds 
to the purely potential velocity field. The exponent < £ < 2 is a free parameter which 
can be viewed as a kind of Holder exponent, which measures "roughness" of the velocity 
field; its "Kolmogorov" value is £ = 4/3, while the "Batchelor" limit £ — > 2 corresponds 
to smooth velocity. The cutoff in the integral ( 12. lip from below at k = m, where 
m = l/£ is the reciprocal of the integral turbulence scale C, provides IR regularization. 
Its precise form is inessential; the sharp cutoff is the simplest choice from calculational 
viewpoints. 

The action functional for the full set of fields $ = ip\ v} becomes 
<S($) = AoVM + {-V t + Ao {d 2 - r ) - a (d iVi )}if) a + 
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It is obtained from (12.71) by the substitution (12 .9p and adding the term corresponding 
to the Gaussian averaging over the velocity field with the correlator (I2.10p . (12. lip : 

S v (y) = -^JdtJd*J dx^x^VH^x'), (2-13) 

where D^ 1 (r) oc D^ 1 r^ 2d ^^ with r = |r| is the kernel of the linear operation inverse to 
Dijir) from (EH]). 

In addition to (12. 8p . the diagrammatic technique for the full-scale model involves 
the velocity propagator (vv) defined by the relations (I2.10p . (12. lip and the new vertex 

^tviV i>ab i/j b = -4>l {(vidi)i/j a + a (diVi)ip a } (2.14) 

with the vertex factor 

V i: a b = -i6 ab (ki + a qi), (2.15) 

where qi is the momentum argument of v j and /c, is the momentum of the field ip. 

As usual for a trilinear interaction, the actual expansion parameter in the model 
(12 .7p is u = gl rather than g itself. Thus for the full model (12. 1 2[) the part of the 
coupling constants (expansion parameters in the ordinary perturbation theory) is played 
by the three parameters 

u = gl ~ A 6 ~ d , w = D /\ ~A^ Wo a ~A § . (2.16) 

The last relations follow from the dimensionality considerations (more precisely, see the 
next section) and set in the characteristic UV momentum scale A. The first relation in 
(I2.16P shows that the interaction tfj^tfj 2 becomes logarithmic (the corresponding coupling 
constant «o becomes dimensionless) for d = 6. Thus for the single-charge problem (12. 7p 
the value d = 6 is the upper critical dimension, and the deviation e = 6 — d plays the role 
of the expansion parameter in the RG approach: the critical dimensions are nontrivial 
for e > and can be calculated in the form of series in e. 

The additional interaction ( 12.14ft of the full model becomes logarithmic for £ = 0. 
The parameter £ is not related to the spatial dimension and can be varied independently. 
For the RG analysis of the full-scale problem it is important that all the interactions 
become logarithmic simultaneously. Otherwise, one of them would be weaker than 
the other from the RG viewpoint and it would be irrelevant in the leading-order IR 
behaviour. As a result, some of the scaling regimes of the full model would be lost. In 
order to study all possible scaling regimes and the crossovers between them, we need a 
genuine multicharge theory, in which all the interactions are treated on equal footing. 
Thus in the following we treat e and £ as small parameters of the same order, e oc £. 
Instead of the ordinary e expansion in the single-charge models, the coordinates of 
the fixed points, critical dimensions and other quantities will be calculated as double 
expansions in e and £. 
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3. Canonical dimensions, UV divergences and renormalization 

It is well known that the analysis of UV divergences is based on the analysis of canonical 
dimensions ("power counting"); see e.g. [H [2]. Dynamical models of the type (12. 7p . 
( 12.12p . in contrast to static ones, have two independent scales: the time scale T and the 
length scale L. Thus the canonical dimension of any quantity F (a field or a parameter) 
is described by two numbers, the momentum dimension d k F and the frequency dimension 
d% defined such that [F] ~ [T]~ d F. By definition 

dk = ~d x = 1, d k = d x = 0, d^ = d t = 0, d UJ = —d t = 1, 

while the other dimensions are found from the requirement that each term of the action 
functional be dimensionless (with respect to the momentum and frequency dimensions 
separately). Then, based on d k F and d F , one can introduce the total canonical dimension 
dp = d F + 2d F (in the free theory, dt oc <9 2 ), which plays in the theory of renormalization 
of dynamical models the same part as the conventional (momentum) dimension does 
in static problems; see e.g. chap. 5 of [2]. The canonical dimensions for the models 
( I2.12p are given in table (TJ including the renormalized parameters (without subscript 
"0"), which will be introduced shortly. 



Table 1. Canonical dimensions of the fields and parameters in the model (|2.12|) . 



F 






V 


Ao, A 


T ,T 


m, /i, A 


u 


w 


u, w, ao, a, a 


4 


(d-2)/2 


(d + 2)/2 


-1 


-2 


2 


1 


6-d 







d F 








1 


1 

















dp 


(d-2)/2 


(d + 2)/2 


1 





2 


1 


6-d 








As already discussed in the end of the preceding section, the full model ( 12.121) 
is logarithmic (all the coupling constants are simultaneously dimensionless) at d — 6 
and £ = 0. Thus the UV divergences in the Green functions manifest themselves as 
singularities in e = 6 — d, £ and, in general, in their linear combinations. 

The total canonical dimension of an arbitrary 1-irreducible Green function T = 
($ ■ ■ • $)i_ir is given by the relation [2] 

d r = 4 + 2dr = d + 2 - N$d$, (3.1) 

where A$ = {N^, N^, N v } are the numbers of corresponding fields entering the function 
T, and the summation over all types of the fields is implied. The total dimension dr in 
the logarithmic theory (that is, at e = £ = 0) is the formal index of the UV divergence 
d~r = dr\ e =(=o- Superficial UV divergences, whose removal requires counterterms, can 
appear only in the functions T for which Sr is a non-negative integer. From table Q] and 
(13~T|) we find 

S T = 8 - 2A^ - 4AV - N v . (3.2) 
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Dimensional analysis should be augmented by certain additional considerations. 
In dynamical models of the type (I2.12p . all the 1-irreducible diagrams without the 
fields vanish, and it is sufficient to consider the functions with Nm > 1 [2]. 
Furthermore, an important role is played by the Galilean symmetry and the invariance 
with respect to the group Q. For example, the function ^ipvv),. can be omitted 
from consideration because the corresponding counterterm iplipaV 2 is forbidden by the 
Galilean invariance. A similar situation occurs for the function {^j wV : the possible 
counterterms ij)\{diVk){dkVj) and ipa(diVk)(diVk) are forbidden by the symmetry with 
respect to Q. In turn, owing to the symmetry with respect to Q, the trilinear term in 
(12.12ft is renormalized as a single entity. 

With those restrictions, the analysis of the expression ( 13. 2 p shows that in our model, 
superficial UV divergences can only be present in the following 1-irreducible functions: 

(^V^i-ir (5 = 0) with the counterterm if) a ip a , 

(tp V)i-ir (8 — 2) with the counterterms ijyjdt^ai Vdc^a, V^a, 

(ijj'ipip) i_ ir (5 = 0) with the counterterm Rabc^l^Pb^c, 
('0t'0i>) 1 _ ir (5 = 1) with the counterterms ip^Vidijipa, tpl(diVi)ip a . 

All these terms are present in the action functional ( I2.12p . so that our model 
appears multiplicatively renormalizable. The Galilean symmetry also requires that the 
counterterms if)'dtif> and ij)*(vd)if) enter the renormalized action only in the form of the 
Lagrangian derivative ^'^ t if), imposing no restriction on the Galilean invariant term 
^(dv)ip. 

We thus conclude that the renormalized action can be written in the form 

S R ($) = ip\ {-Z.Vt + A (Z 2 d 2 - Z 3 r) - aZ 6 (d iVi )}i> a + 

+ XZ^l - g^l 2 R abc \Z^J b i; c /2 + S v (v) (3.3) 

with the same S v (v) from ( 12. 131) . 

Here A, r, g and a are renormalized counterparts of the bare parameters (with the 
subscripts "0" ) and \i is the reference momentum scale (additional arbitrary parameter 
of the renormalized theory). The renormalization constants Z absorb the singularities 
in e and £ and depend on the dimensionless parameters u, w, a and a. Expression 
( 13. 3p can be reproduced by the multiplicative renormalization of the fields ip — > ipZ^, 
—> ip^Z^ and the parameters: 

£o = gV E/2 Z g , u = ufi e Z u , w = WfjtZy,, 

A = \Z X , r = tZ t , a = aZ T . (3.4) 

Since the last term S v (y) in ( 13. 3 p remains intact, the amplitude Do from (12.111) is 
expressed in renormalized parameters as D = w Xo — w\fit, while the parameters m 
and a are not renormalized: mo = m, «o = ct. Owing to the Galilean symmetry, the 
both terms in the covariant derivative Vt are renormalized with the same constant Z\, 
so that the velocity field is not renormalized, either. Hence the exact relations 

Z W Z\ = 1, Z m = Z a — Z v — 1. (3-5) 
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{^V)l-ir = 2XZ S + - 





{^h-ir = i^^i - a izy + z 3T 





A A 




+ 




•♦A- A -A- A 



Figure 1. Superficially divergent 1-irreducible Green functions in the one- loop 
approximation. 



Comparison of the expressions (12.1 2p and (13.31) gives the following relations between the 
renormalization constants Z\-Z§ and (13. 4|) : 



Z x 
Z* 



Z%pZ^, z 2 

Z\ZgZ^Z^ , 



ZiZ, 



ZxZ 



Z 2 Z T , 

Zq = Z\Z a . 



(3.6) 



Resolving these relations with respect to the renormalization constants of the fields and 
parameters gives 



Z\ Zi, Zf — 

7 l/2 7 l/2 7 -l/2 
^1 ^2 ^5 > 



Z 2 1 Z3, 



Z a — Z l Zq, Z u — Z x Z 2 Z^Z^, 

7 l/2 7 -l/2 7 l/2 
^1 ^2 ^5 > 



(3.7) 



7 2 

Z 9- 



^ ~~ ^1 ^2 ^5 ' Z; V' t ~~ ~1 ^2 

where we have passed to the coupling constant u = g 2 with Z u 

The renormalization constants can be found from the requirement that the Green 
functions of the renormalized model ( \3.3h . when expressed in renormalized variables, be 
UV finite (in our case, finite at e — > 0, £ — > 0). The constants Z\-Z§ are calculated 
directly from the diagrams, then the constants in (13. 4p are found from (13. 7p . In order to 
find the full set of constants, it is sufficient to consider the 1-irreducible Green functions 
which involve superficial divergences. The diagrammatic representation for the relevant 
Green functions in the one-loop approximation is given on figure [U 

The solid lines with arrows denote the propagator (■0^') o , the arrow pointing to 
the field The solid lines without arrows correspond to the propagator (ipip)o and 
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the wavy lines denote the velocity propagator (vv) defined in (I2.10p and (12. lip . The 
external ends with incoming arrows correspond to the fields ifi\ the ends without arrows 
correspond to ip. The triple vertices with one wavy line correspond to the vertex factor 
(I2.15p . In order to calculate the renormalization constants, it is sufficient to consider 
the symmetric phase with r > 0. Then, with respect to the "isovector" indices, all the 
terms in the first three lines on figure 1 are proportional to 5 a i, while all the terms in 
last line are proportional to R a bc- 

All the diagrammatic elements should be expressed in renormalized variables using 
the relations ( I3.4p -( l3~6j) . In the one-loop approximation, the Z's in the bare terms 
should be taken in the first order in u = g 2 and w, while in the diagrams they should 
simply be replaced with unities, Z — > 1. Thus the passage to renormalized variables in 
the diagrams is achieved by the simple substitutions Ao — > A, To — » r, go — > g/i £//2 and 

Wo — > 

In practical calculations, we used the minimal subtraction (MS) scheme, where 
all the renormalization constants have the forms U Z = 1+ only singularities in e 
and with the coefficients depending on the completely dimensionless renormalized 
parameters u, w, a and a. 

The one-loop calculations for similar models are discussed in detail, e.g., in [T9]-|21|. 
and here we only give the results: 

Zi = l-——, Z 2 = l-— — (5 + a), Z 3 = l , 

2e 3e 64 e 

n 1 2R 2 u waa 2 uRi w . lN2 

Z 4 = 1 — , Z 5 = 1 - — — a(a - 1) , 

£ £ 2e £ 

Z 6 = l-^i(4a-l) (3.8) 
Zae 

with Rip from ( 12.4ft and up to the corrections of the order u 2 , w 2 , uw and higher. To 
simplify the resulting expressions, we have passed to the new parameters 

u -)> u/1287r 3 , w -> w/647r 3 ; 

in (13.8p and below they are denoted by the same symbols u and w. 



4. RG functions and equations 

Let us briefly recall an elementary derivation of the RG equations; detailed presentation 
can be found e.g. in [TJ [2]. The RG equations are written for the renormalized 
Green functions W R = (<&••■<!>)#, which differ from the original (unrenormalized) ones 
W = only by normalization (due to rescaling of the fields) and choice of 

parameters, and therefore can equally be used for analyzing the critical behaviour. The 
relation Sr(Q, e, /x) = S(Z$$, cq) between the bare ( I2.12p and renormalized ( 13.31) action 
functionals results in the relations 

W(e ,...)=Z^Z*fW R (e, »,...). (4.1) 
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between the Green functions. Here and are the numbers of corresponding fields 
in the function W (we recall that in our model Z v = 1); e = {A , r , u , w , a , m , a } 
is the full set of bare parameters and e = {A, r, u, w, a, m, a} are their renormalized 
analogs (we recall that = a and m = m); the ellipsis stands for the other arguments 
(coordinates/momenta or times/frequencies). 

We use £> M to denote the differential operation fid^ for fixed e and operate on both 
sides of the relation (14. ip with it. This gives the basic RG differential equation: 

{V RG + + iV>7^t } W R (e, fi,...) = 0, (4.2) 

where T>rg is the operation D M expressed in the renormalized variables: 

V RG = V^ + Mu + + (5 a d a - 7A £> A - lr V T (4.3) 

and we have written T> x = xd x for any variable x. The anomalous dimensions 7 are 
defined as 

7i? = In Zp for any quantity F, (4.4) 

while the /3 functions for the coupling constants u, w and a are 

U = V^u = u{-s- 7 U ), 
(3 W = V^w = w (-£ - j w ), 

a = V^a = -aj a , (4.5) 

where the second equalities come from the definitions and the relations (13 .4ft . In 
principle, the dimensionless parameter a should be treated as the fourth coupling 
constant, but the corresponding (3 function 

(3 a = V^a = -a^f a (4.6) 

vanishes identically due to (I4.12p and does not appear in the subsequent relations. 

The anomalous dimension corresponding to a given renormalization constant Zp is 
found from the relation 

If = (AA + P w d w + (5 a da) In Z F ~ - (eV u + £V W ) In Z F . (4.7) 

The first equality follows from the definition (I4.4p . expression (14. 3p for the operation 
P M in renormalized variables, and the fact that the Z's depend only on the completely 
dimensionless coupling constants u, w and a. In the second (approximate) equality, we 
only retained the leading-order terms in the functions ( I4.5j) . which is sufficient for the 
first-order approximation. The leading-order expressions (13 .8p for the renormalization 
constants have the forms 

11 w 

Z F = 1 + -A F (a,a) + —Bjr(a,a). (4.8) 

The factors e and £ in (14 .7\\ cancel the corresponding poles contained in the expressions 
(14. 8 p for the constants Z F , which leads to the final UV finite expressions for the 
anomalous dimensions: 

j F = — uA F (a, a) — wB F (a, a) (4.9) 
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for any constant Z F . Then equations ( 13. 8 p give 

7i = R x u/2, 72 = -Riw/3 + w(5 + a)/6, 73 = 2_Riw, 
74 = 2i? 2 ^ + w«a 2 , 75 = uR\/2 + wa(a — l) 2 , 

7e = ui2i(4a-l)/2a. (4.10) 

The multiplicative relations (13. 6p between the renormalization constants result in 
the linear relations between the corresponding anomalous dimensions: 

7a = 72 - 7l> It = 73 - 72, 

7a = 76 - 7i, 7« = -7i - 372 + 274 + 75, 

27v = 7i + 72 - 75, 27^=71-72 + 75, (4.11) 
while the exact relations (13.51) result in 



lw = -7A, 7m = la = lv = 0. (4.12) 

Along with (14.101) . those relations give the final first-order explicit expressions for the 
anomalous dimensions of the fields and parameters: 

7a = — j w = — Riu/6 + w(5 + a)/6, 7 T = 5Riu/3 — w(5 + a)/6, 
7 M = — uR — w(B + a)/2 + watf(a), 7 a = (3a — l)Riu/2a, 
27^, = Rxu/3 — w [5w(a — l)/6 + «a(a — 2)] , 

2 7v; t = 2i? 1 u/3 + [5w(a - l)/6 + aa(a - 2)] , (4.13) 
where we have denoted R = R\ — 4i?2 and f(a) = 2a 2 + (a — l) 2 . 

5. Attractors of the RG equations and scaling regimes 

It is well known that possible asymptotic regimes of a renormalizable field theoretic 
model are determined by the asymptotic behaviour of the system of ordinary differential 
equations for the so-called invariant (running) coupling constants 

U s gi{s,g) = pi(jg), gi(l,g) = g h (5.1) 

where s — k/fj,, g — {gi} is the full set of couplings and (ji(s,g) are the corresponding 
invariant variables. As a rule, the IR (s — » 0) and UV (s — > 00) behaviour of such 
system is determined by fixed points g^. The coordinates of possible fixed points are 
found from the requirement that all the (3 functions vanish: 

AG?*) = 0. (5.2) 

The type of a given fixed point is determined by the matrix 

= <9/3; /<%!<?=<?* , (5.3) 

which appears in the linearized version of the system ( 15.11) near that point. For IR 
attractive fixed points (which we are interested in here) the matrix Q is positive, i.e. 
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the real parts of all its eigenvalues are positive. In the case at hand, the fixed points for 
the full set of couplings u, w, a, a are determined by the equations 

with the functions defined in the previous section. However, in our model the 
attractors of the system (15. ip involve, in general, two-dimensional surfaces in the full 
four-dimensional space of couplings. Indeed, the function f 14 . 6 1) vanishes identically, and 
the equation (3 a = imposes no restriction on the parameter a. Thus it is convenient 
to consider the attractors of the system ( 15.11) in the three-dimensional space u, w, a; 
their coordinates, matrix (15. 3p and the critical exponents depend, in general, on the free 
parameter a. Furthermore, in this reduced space some attractors will be not simply 
fixed points, but also lines of fixed points, parametrized by the coupling a. In the 
following we will use the term "fixed point" for all those attractors, bearing in mind 
that their coordinates can depend on a and, in general, on a. 

The couplings and w* should be positive (by definition, u oc g 2 > and 
w oc Dq/X > 0), so that the point is admissible from the physics viewpoints if it 
satisfies the conditions 

> 0, > (5.5) 

and can be IR attractive (Q > 0) for some values of the model parameters (in fact, all 
the above inequalities can be non-strict). 

In the one-loop approximation, the /3 functions are found from the definitions ( 14.5P 
and the explicit expressions ( 14.121) and ( 14. 13ft : 

(3 U = u [—e + Ru + w(5 + a)/2 — waf(a)} , 
W — w [— £ — Riu/6 + w(5 + a)/6] , 

a =uR 1 (l-3a)/2, (5.6) 

with R — R 1 — AR 2 and f(a) = 2a 2 + (a - l) 2 . 

The general pattern of the attractors of the system ( 15.61) is rather complicated and 
depends qualitatively on the values of the parameters i?i j2 and a. There are four fixed 
points: 

(I) The line of Gaussian (free) fixed points: = = 0, a* arbitrary. 

(II) The point with = 0, corresponding to the pure Potts model (turbulent 
advection is irrelevant). 

(III) The line of fixed points with u* = and arbitrary a* corresponding to the 
passively advected scalar without self-interaction. 

(IV) The most nontrivial fixed point with 7^ 0, w* 7^ 0, corresponding to the 
novel scaling regime (universality class): both the advection and the self- interaction are 
relevant. 

Let us discuss these points in more detail. 
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5.1. Fixed points with = 0. 

For the Gaussian fixed point (point I) u* = = 0, a* arbitrary. The only nonzero 
off-diagonal element in Q is Q au . Thus the matrix Q is triangular and its eigenvalues 
coincide with the diagonal elements: Q u = —e, Q w = — £, Q a = 0. Here and below 
Qi = Qa denote the diagonal elements (no summation over i). Vanishing of the last 
eigenvalue reflects the fact that the point is degenerate. 

The coordinates of the passive scalar point (point III) are: 

xl* = 0, = r, a* arbitrary. (5.7) 

(5 + a) 

Now Q ua = Q wa = Q uw = 0, so the matrix Q is block-triangular, and the eigenvalues are 
given by the diagonal elements: 

Q a = 0, Q W = Z, n u = -e + 3Z-f(a)j^- (5.8) 

(5 + a) 

The condition Q w > also gives > 0. 

The function f(a) = 2a 2 + (a — l) 2 achieves the minimum value /(1/3) = 2/3 at 
a = 1/3, and in Q u we can write f(a) = /(1/3) + If (a) — /(1/3)]. This gives 

V u = 0^=1/3 - 7^r[/(a) - /(1/3)] > 0. (5.9) 
(5 + a) 

The second term is negative, so (15. 9p can be satisfied only if f2J a =i/3 > 0. This gives 

- (5 + a)e + (15-a)^ > 0. (5.10) 

This is the domain in the £— e plane where the point can be stable. Now ( 15. 9ft gives the 
restriction for a*: 

6<[/(«) - /(1/3)] < -(5 + a)e + (15 - a% (5.11) 

which gives 

( i ^ -(5 + a)e + (15-a)^ 

(a* - 1/3 < — — . 5.12 

18a4 

We conclude that the admissible fixed points of the type III form an interval on the 
line ( 15 .71) specified by the inequality (15.121) . The region of IR stability in the £-e plane is 
given by the inequalities ( I5.10p and ^ > 0, so that the condition w* > is automatically 
satisfied. 

For a = the boundary defined by the inequality (15.101) is £ = e/3. When a 
increases, it rotates counter clockwise, and for a — > oo tends to £ = — e. 



5.2. Fzxec? points with u* ^ 0. 

Now let us turn to the fixed points with u ^ 0. Then the equation /3 a = readily gives 
a = 1/3. Furthermore, one has f2 ua = Q wa = 0, and the eigenvalue Q a = —SRiu/2 
decouples. For u > 0, it can be positive only if R\ < 0. Thus in the following we assume 
Ri < 0; the case R\ > requires special attention and will be discussed later. 
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We put a = 1/3 in ft u>w and arrive at a closed system of two ft functions for the 
two couplings u, w of the form: 

ft u = u[-e + Au + Bw], ft w = w[-£ + Cu + Dw}. (5.13) 

For our set of ft functions (15.61) the coefficients in ( 15 . 13f) are: 

A = R, B = (15 - a)/6, C = -Ri/6 > 0, D = (5 + a)/6 > 0, (5.14) 

but it is instructive to discuss it first in a general form with arbitrary real coefficients 
A-D. Now we are interested only in the fixed points with « ^ 0; there are two such 
points: the Potts-type point with w = and the full-scale point with 

Pure Potts point (point II). Here u>* = 0, = e/R. Now Q wu = and the reduced 
2x2 matrix Q is triangular. Then the point is IR stable for e > (so that R > 
since we require that u* > 0) and — £ + > 0. The last relation gives £ < —eR\/6R, 
because A/C > 0. Thus point II 

= e/H, u>* = 0, a* = 1/3 

can be physical only if R > 0, Ri < (and any a) and is IR stable in the region 

e > 0, £ < —Rie/QR. (5.15) 

The coordinates of the full-scale fixed point (point IV) are 

u* = (De-Bg)/A, w* = (A£-Ce)/A 1 A = AD — BC, (5.16) 

while the reduced matrix Q can be written in the form 

It is useful not to substitute explicit expressions (15.161) for a while. For the positive 2x2 
matrix Q the eigenvalues can be real and positive, or they can be complex conjugate 
with positive real parts. Thus the necessary and sufficient condition for the IR stability 
can be is given by the two inequalities: 

dettt > 0, trft > 0. (5.18) 

From (I5.17P we obtain 

det Q = u*w* A > 0, (5.19) 

which along with ( 15. 5 j) shows that this point can be admissible only if A > 0. For trfl 
we obtain: 

trfi = Au, + Dw* > 0. (5.20) 

There are three possibilities: 

(1) A > 0, D > 0. In this case the inequality ( 15.20P is an automatic corollary of 

flESD. 

(2) A < 0, D < 0. Then f l5.20f) contradicts to ( 15. 5ft and this point cannot admissible. 

(3) The parameters A and D are opposite in sign. For definiteness, we assume that 
A < 0, D > 0. In this case from ( I5.20p one obtains: 

w* > -AujD > 0, (5.21) 
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where the last inequality follows from A < and > 0. The second inequality in ( 15. 5p 
is implied by the ( I5.2ip and thus becomes superfluous. The region where the fixed point 
is IR attractive and positive is given by the two inequalities 

> 0, Au* + Dw* > 0. (5.22) 

We conclude that point IV can be admissible if A > and at least one of the 
two parameters A and D is positive. The region of admissibility is determined by the 
inequalities ( 15. 5ft if A and D are both positive, and by the inequalities of the type (15. 22ft 
if A and D have different signs. 



Now we turn to our specific model with the j3 functions (15. 6p and the coefficients 
( I5.14p . Then the coordinates of the fixed point IV have the forms 

Ms + oO-jtis-aQi W 6 + Ja a- = 1A (5 . 23) 

oA A 

where 

5 ct 

A = AD — BC = -(R + Ri 2) + -(R - R 1 6 . (5.24) 
o o 

This point can be physical only if A > 0. Since D > 0, A = R can be of either 
sign; we also recall that R± < 0. There are four different cases: 

(z) R > -Ri/2, 

(ii) - Rx/2 > R > 0, 

(Hi) > R > i?i/6, 

(iv) R < Rt/6, (5.25) 
which should be discussed separately. 



Case («). Then automatically R > and R — R\/6 > 0, so that A > for all a. 
The conditions that point IV is IR attractive coincide with the conditions ( 15. 5ft that its 
coordinates are positive. 

Case (ii). Then R — i?i/6 > 0. Thus A < for small a, but becomes positive for 
a > a , where 

= (R + Ri/2) 

(R-Rt/6) v ; 

The conditions that the point is IR attractive are again ( 15.51) . 

Case (Hi). Then R — i?i/6 > 0, and the point IV can be admissible for a > a with 
the same ao- Now A < 0, and the conditions that the point is physical are given by the 
inequalities (15.221) . 



Case (iv). Then A < for all a, and point IV cannot be admissible. 
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Now let us write the positivity conditions ( 15. 5 p for the cases (i) and (ii) in a more 
explicit form with the aid of expressions (15.161) and (15.231) : 

{A£ - Ce) > 0, (De - B£) > 0. (5.27) 

Since R = A > 0, the first inequality is £ > (C/A)e. Since B > for a < 15, the second 
inequality is £ < (D/B)e. It is also important that 
C D —A 



A B AB 

so that 



< (5.28) 



C/A < D/B. 

Also note that C/A and D/B are positive. 

For a > 15 we have B < 0, and the second inequality becomes £ > (D / B)e with 
D/B < 0. 

Thus for a < 15, the admissibility region is the sector in the upper right quadrant 
in the e— £ plane, bounded by the ray £ = (C/A)e from below and £ = (D/B)e from 
above. When a grows, the upper ray £ = (D/B)e rotates counter clockwise and moves 
to the upper- left quadrant. For case (i) a changes from to oo and the ray changes 
from £ = e/3to£ = — e (exactly like the boundary (I5.10p of point III). 

For case (ii) a changes from a to oo and the ray changes from £ = —eR\/QR to 
£ = — e. For the following it is important that for case (ii) one has —R\/6R > 1/3. Also 
note that at a = «o the two boundaries for point IV coincide with each other (this is 
also obvious from relation ( 15. 28ft . in which A = for a = cuo) and with the boundary 
(I5.15P of point II. Also note that 

a °- 15 =(S^6) <0 - (5 ' 29) 

Let us turn to the case (Hi). Now the inequality u* > in (I5.22p becomes 
B£ < De. 

Now from ( I5.29P we see that «o — 15 > 0, so that B < and we obtain 

i > e(D/B) with D/B < 0. (5.30) 

The second condition Au* + Dto» > in ( I5.22p is 

iA(D — B) > sD(C - A) 

where A = R < and (D — B) = (a — 5)/3 > 0. The last inequality holds because 
a > cto > 5: 

(R + RJ&) 
a °- 5 = - 10 (R-R 1 /6) > °- 
Thus the second inequality is 
D(C-A) 

« < e io^wr (5 - 31) 
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where 

D{C -A) (5 + a)(R + R x /$) 



< 0. (5.32) 



A(D - B) 2R(5 - a) 

It is also important that D/B > D{C — A)/ A{D — B), because 
D(C-A) D —AD ^ 

A(D -B) B~ AB(D - B) < U ' ^ 66} 

Thus the both inequalities ( I5.30f> . (I5.3ip are satisfied in a sector in the upper left 
quadrant; the lower bound is f)5.30p and the upper bound is (I5.3ip . 

From ( 15.321) it follows that, when a changes from ao to oo, the coefficient 
D(C - A)/A(D - B) changes from -RJ6R to -(R + R 1 /6)/2R. The coefficient 
D/B = (5 + a)/ (15 — a) changes from —Ri/QR to —1. Thus for a = a® the domain has 
zero width (this is also obvious from the relation (15. 33ft . in which A = for a = ceo), 
and when a increases it is getting wider. 



5.3. General pattern of the fixed points. 

Now we are in a position to describe the general pattern of the admissibility regions of 
the fixed points in the e-£ plane. In the one-loop approximations, they all are sectors 
bounded by straight rays; in the following, they are referred to as sectors I, II etc. There 
are four different situations related to the four cases in (I5.25p . 

Case (i) is illustrated by fig. [2j The e-£ plane is divided into four sectors LTV 
without gaps or overlaps. The boundary between sectors II and IV (solid line) is given 
by the ray £ = —sRi/6R. For a small, the boundary between sectors III and IV (dashed 
line) lies in the upper right quadrant (fig. |2]a). As a grows, it rotates counter clockwise 
and for a > cto moves to the upper left quadrant (fig. E&). Here and below, dotted lines 
denote the a — > oo limits of various boundaries. The boundary between II and IV (solid 
line) is given by the ray £ = —eRi/QR. 

Case (ii) is illustrated by fig. [3j For a small, sector IV is absent, while sectors I III 
cover the entire plane e-£ without gaps, but with an overlap between II and III (fig. [3]a): 
the boundary £ = —eRi/6R of sector II (solid line) lies above the boundary £ = e/3 
of sector III (dashed line). The existence of overlap means that, for the corresponding 
values of e and £, the critical behaviour is non- universal: it can be described by the fixed 
points II or III, depending on the initial data for the problem (15.11) . As a grows, the 
boundary of sector III rotates, the overlap is getting thinner and disappears for a = olq. 
For a > cto, there is a gap between sectors II and III. Meanwhile, sector IV appears and 
it fills that gap exactly, (fig. E]6) . 

Thus starting with a = a , the £-£ plane is divided into four sectors LTV with no 
gaps nor overlaps. 
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Figure 2. Case (i); a — (left) and a > ao (right). 




Figure 3. Case (ii); a — (left) and a > ao (right). 



Case (Hi). Sector II is absent, while IV is absent for small a. There is an empty 
space in the e— £ plane, not covered by any of the admissibility sectors (fig. H]a). This is 
interpreted as absence of a second-order phase transition for such values of e and £. As a 
grows, the boundary of sector III (dashed line) rotates counter clockwise, and at a = a 
sector IV, adjacent to III, appears in the upper left quadrant. Its right boundary (solid 
line) also rotates counter clockwise, but such that its width increases with a. There is 
no gap between III and IV for all a > ao (fig. H]6). 

Thus the turbulent mixing can lead to the emergence of a critical state in a situation, 
where admissible fixed point does not exist for the original static model (I2.3P and for 
the equilibrium stochastic problem (12. lft . (12. 2 ft without mixing. 



Case (iv). The most "boring" case, illustrated by fig. El Sectors II and IV are 
absent for all a. When a grows, sector III decreases, while the empty region grows. 
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Figure 4. Case (Hi): a = (left) and a > ao (right). 




Figure 5. Case (iv). 



5.4- Fixed points with R\ > 0. 

It remains to discuss the case R\ > and u* ^ 0. Then the fixed points with a* = 1/3 
cannot be IR attractive; see the discussion in section 15.21 However, nontrivial points for 
that case can be found in terms of the new couplings b = 1 /a and v = wa 2 with the /3 
functions 

A, = (-l/a 2 )(3 a , (3 V = a 2 (3 w + 2aw(3 a . (5.34) 

Then 

p b = -uRtbib - 3)/2. (5.35) 

The relevant fixed point is 6* = with = Qbu = 0, so that fi& = 3uRx/2 > (for 
u, R\ > 0) is an eigenvalue, and b decouples. We put b = in the other (3 functions and 
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again obtain a closed system of the type (15.131) : 

p u = u[-e + Ru - 3av], f3 v = v[-£ - 19ui2i/6]. (5.36) 

We immediately see that A = — (19/2) ai?i < 0, so that the full-scale point with u^O, 
w 7^ cannot be admissible. The other possible fixed point is 

f* = 0, it* = e/R, (5.37) 

with the IR stability conditions 

e > 0, f < -19i2ie/6i2, (5.38) 

so that R > 0. Then the sector of admissibility lies in the lower right quadrant. The 
point is similar to II in the sense that v * = and the advection is irrelevant. 

5.5. The plain Potts model 

Let us conclude this section with a brief discussion of the original Potts model with 
the hypertetrahedron symmetry. From (12.61) we obtain R = (n + 1) 2 (7 — 3n). The 
most interesting cases are n = (percolation process in a moving medium) and n = 2 
(nematic-to- isotropic transition in a liquid crystal). For n = we have Ri = — 1, 
R 2 = —2 and R = 7; those values belong to case (i) from (15.251) . For incompressible or 
weakly compressible fluid, the most realistic values £ = 4/3 and g? = 3 correspond to 
the passive scalar regime (point III). As a increases, the boundary between the regions 
III and IV moves, and for a large enough the same values correspond to the new regime 
(point IV). Thus the compressibility leads to the changeover in the type of critical 
behaviour between two universality classes. 

For n = 2 we have R 2 = and R\ = R = 9. Thus the full-scale fixed point 
cannot be admissible, while the regions of admissibility of the Potts-type point II and 
the passive scalar case with 6* = lie in the lower right quadrant (and thus are not very 
interesting for physics applications). For small a, the aforementioned physical values 
of £ and e belong to the passive scalar case with finite a* (point III), while for a large 
enough they correspond to neither admissible point. Here the growth of compressibility 
destroys the critical state. 

For n > 3 we have R\ > and R < 0, so that the points with 6* = and the 
Potts-type point II cannot be admissible. Depending on the value of a, the physical 
values of £ and e belong to the passive scalar case ( 15.371) or lie in the "desert" in the 
£-£ plane with no admissible fixed points. 

6. Critical scaling and critical dimensions 

Existence of IR attractors in the RG equations implies existence of asymptotic scaling 
regimes for all the Green functions in the IR range. In dynamical models, critical 
dimensions Ap of the IR relevant quantities (times/frequencies, coordinates/momenta, 
t and the fields) F are given by the relations (see e.g. chap. 5 in [2]) 

A F = d k F + A u d F + lF , A w = 2- 7 *. (6.1) 
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Here d!p U are the canonical dimensions of F, given in table [I], and 7^ is the value of 
the corresponding anomalous dimension ( 14. 4 p at the given fixed point. In our case, 
If — 1f(u*, w*, a*, a). This gives: 

A^ = (d-2)/2 + %, A^ t = (d + 2)/2 + 7 ; t , A r = 2 + 7 ;. (6.2) 

The final results are obtained by substituting the coordinates of the fixed points into 
the explicit one-loop expressions (14. 13j) for the anomalous dimensions. 

In particular, the response (Green) function in the IR range takes on the asymptotic 
form (in the symmetric phase r > 0) 

(^ a (t,x)Vi(0,0)> = 6 ab r' A ^ A ^ $ (tr- A ^rr^) (6.3) 

with r = |x| and some scaling function $. 

For the Gaussian fixed point I the dimensions are trivial: 7^ = for all F (here and 
below, we present the results for the anomalous dimensions, which are more graspable). 

For point II the known one-loop results for the Potts model are recovered: 

* _ Ri£ * Ri£ * _ 5_Ri£ *_ _^i£ (f- a\ 

7v» - 6R > 7^t - 3R , 7 r ~ 3R , 7a - QR , 

with corrections of order 0(e 2 ) and higher (no dependence on £). 
For the passive scalar point III one derives: 

(5 -5a- Qaal + 12aa*) , 
7^ = -7,t = C, 7A = -7 T = e, (6-5) 

where a* lies in the interval (I5.12p . The expressions for 7 ^ r , as well as the relation 
7^ = — 7^ t (so that A^ + A^t = d), are exact, because they actually refer to the usual 
Kraichnan's model without self-interaction. 
For the full-scale point IV we obtain: 

7 * t = (45i?ie + 17.Ri as - 180i2i£ + 12i2ia£ - 90i^ + 30ito£)/216A , 

7 ; = (45 J R 1 £ + P^ae - 90R^ + 6Ria£ + 9QR£ - 30itof)/216A , 

7* = (45i?ie + 9i?iae - 150i?i£ + 10i?i< - 30i?£ - 6i?c<)/36A (6.6) 

with A from ( I5.24p and the higher-order corrections in £ and e. The exact result 7^ = £ 
follows from the first relation in (14. 12|) and the equation fi w = with w* ^ 0. 



7. Conclusion 



We studied effects of turbulent mixing on the critical dynamics of a nearly critical 
system, whose equilibrium behaviour is described by the Ashkin- Teller-Potts model. 
The turbulent mixing was modelled by Kraichnan's rapid-change ensemble: time- 
decorrelated Gaussian velocity field with the power-like spectrum cx k~ d ~^. Special 
attention was paid to compressibility of the fluid, because it leads to interesting 
qualitative crossover phenomena. 

The original stochastic problem was reformulated as a multiplicatively renormal- 
izable field theoretic model, which allowed us to apply the field theoretic RG to the 
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analysis of its IR behaviour. We showed that, depending on the relation between the 
space dimension d, the exponent £ and the degree of compressibility, the model reveals 
four types of possible IR asymptotic behaviour, associated with the four attractors (fixed 
points) of the RG equations. Three fixed points correspond to known regimes: Gaus- 
sian (free) theory, passively advected scalar field and the original Potts model without 
mixing. The most interesting fixed point corresponds to a new type of critical behaviour 
(universality class), where the self-interaction of the order parameter and the turbulent 
mixing are equally important, and the critical dimensions depend on d, £, the symmetry 
group Q and the compressibility parameter a. 

Explicit results were derived within the leading (one-loop) approximation, that is, 
in the leading order of the double expansion in e and £. Thus their validity for finite 
physical values of these parameters can be called in question, especially because of 
large physical values e = 6 — d = 3 and £ = 4/3. Careful discussion of this problem 
requires analysis of higher-order corrections and applying some kind of resummation 
procedure. Such analysis goes far beyond the scope of the present paper, and we hope 
to address it in the future. Nevertheless, the discussion of the RG flows, given in [22] 
for a similar problem, suggests that the pattern of critical regimes, obtained in the one- 
loop approximation, appears robust with respect to higher-order corrections and can be 
preserved for finite values of e and £. 

Further investigation should account for conservation of the order parameter, its 
feedback on the dynamics of the velocity statistics, finite correlation time and non- 
Gaussian character of the advecting velocity field. This work is in progress. 
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